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ABSTRACT 

We report the results from the temporal and spectral analysis of an XMM-Newton 
observation of Nova Centauri 1986 (V842 Cen). We detect a period at 3.51±0.4 h in the 
EPIC data and at 4.0±0.8 h in the OM data. The X-ray spectrum is consistent with 
the emission from an absorbed thin thermal plasma with a temperature distribution 
given by an isobaric cooling flow. The maximum temperature of the cooling flow model 
is kTmax — 43l]^2 keV. Such a high temperature can be reached in a shocked region 
and, given the periodicity detected, most likely arises in a magnetically-channelled 
accretion flow characteristic of intermediate polars. The pulsed fraction of the 3.51 h 
modulation decreases with energy as observed in the X-ray light curves of magnetic 
CVs, possibly due either to occultation of the accretion column by the white dwarf 
body or phase-dependent to absorption. We do not find the 57 s whi t e dwa rf spin 
period, with a pulse amplitude of 4 mmag, reported by IWoudt et al.l ()2009D either 
in the Optical Monitor (OM) data, which are sensitive to pulse amplitudes > 0.03 
magnitudes, or the EPIC data, sensitive to pulse fractions p >14±2%. 
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1 INTRODUCTION 

Cataclysmic variables (CVs) are binary systems consisting of 
a white dwarf primary and a secondary that transfers mass 
to it through an accretion disk around the white dwarf. Their 
brightness can change by large factors (several millions) dur- 
ing nova outbursts, which are due to thermonuclear run- 
aways of the hydrogen-rich material that has accreted onto 
the white dwarf. 

The power spectra of CVs display several different pe- 
riodicities which are attributed to the orbital period, white 
dwarf spin period (in the case of magnetic CVs) and su- 
perhumps (which are photometric modulations near, but 
not at, the orbital period) among others. Among magnetic 
CVs, intermediate polars (IPs) are moderately strong X-ray 
sources {Lx ~10'^^~'^'^ erg s~^), in which accreting mate- 
rial is magnetically channeled and shock-heated in the ac- 
cretion column, producing high energy radiation as it cools 
before reaching the white dwarf surface. Strong X-ray mod- 
ulation is observed at the spin period, likely due to self- 
occultation of the accretion column by the white dwarf's 
body and/ or phase-dependen t absorption by pre-shock ma- 
terial (e.g.. I Allan et al.lll998h . 

Nova Centauri 1986 (V842 Cen) underwent its lat- 
est nova outburst in February 1986 and took 48 days 



to decline its brightness by 3 magn itudes, which suggest s 
that it was a m oderately fast nova 'S ekiguchi et al.lll989l ). 
iDownes fc Due rbcck (2000) imaged the nova ejecta in 1998 
using H/3 and [O IIIJA5007 narrowband filters. The ejecta 
show a fairly spherical shape in both filters (5.6" x 6.0"). 
Given this information, therefore, nothing is particularly pe- 
culiar about this nova. However, optical photometry obser- 
vations in 2008 overturned this picture when a fast p ulsation 
of 57 seconds was discovered bv lWoudt et al.l (|2009l ') and at- 
tributed to the spin period of the white dwarf, which would 
make V842 Cen a member of the IP class. This finding has 
strong implications: its white dwarf would be the fastest ro- 
tator currently known in a nova remnant; and, its symmetric 
shell morphology would now present a puzzle, since a mag- 
netic field anchored to a rapidly rotating white dwarf would 
be like ly to produce asymme t ric ejecta sh apes during a nova 
event (JFiedler fc Jonedll980l : lLiviolll995l '). 

In this letter we show the results of an 'K}A}A-N ewton 
X-ray and optical observ ation aimed to sear ch for the 57 s 
spin period discovered bv lWoudt et al.l ((20091) and thus con- 
firm the IP nature of V842 Cen. In Section [2] we present 
the details of the observation and data processing, while in 
Section [3] we show the results obtained from the timing and 
spectral analysis. Discussion and conclusions are presented 
in Section |4l 
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2 OBSERVATIONS AND ANALYSIS 

We observed V842 Cen with the IQAyi- Newton Observatory 
on 2011 February 24 for 56.9 ks using the EPIC instrument, 
operated in full window mode with the medium thickness 
filter and the Optical Monitor (OM) in fast mode. After 
removing events at periods with high flaring particle back- 
ground using a Scr clipping method, the resulting exposure 
time reduced to 49.8 ks. 

The source spectra and light curves were accumulated 
from circular regions of 32" and 20", respectively, centered 
on V842 Cen (a=14h 35m 52.55s, (S=-57° 37' 35.3"). The 
background spectra and light curves were extracted from a 
source-free region on the same chip taken within a circle of 
40" radius. The total number of counts in the source region 
in the pn, MOSl and M0S2 cameras was 5354 which in- 
cludes 735 counts from background. We used the rmfgen 
and arfgen to build the redistribution matrices and an- 
cillary responses, respectively. The resulting X-ray spectra 
were grouped with a minimum of 15 counts per energy bin. 
For timing analysis, we converted photon arrival times to 
the solar system barycenter using the SAS task barycen. 
The exposure times in the OM filters were 16594, 16599 and 
16596 seconds for the V, U and B respectively. The OM light 
curves were reprocessed using the script omf chain with time 
bin size of 0.5 s. 



3 RESULTS 

3.1 Timing 

For the X-ray data, we used the event arrival times and 
Rayleigh statistics (Zf; iBuccheril [l983 : IStute et al] I2OIII I 
to search for periods in the frequency range from fmin = 



^ / J- span rO Jrr 



l/2tfr 



with A/ = 1/ATspan, where 



1^ frame IS me 

EPIC camera readout time (73 ms) and A is the oversam- 
pling factor, which we set equal to 1,000. We detect a period 
at 3.64±0.46 h which is within the errors of the ~3.94 h or- 
bital period (Fig. [T]) estimated bv lWoudt et all (|2009l ). We 
are not able to detect the ^^57 s white dwarf spin period in 
the X-ray data. 

Since the particle background is not expected to be 
modulated, it should be safe to use all the exposure time 
(56.9 ks) to improve the accuracy of the period detection. 
In this case, the modulation is found to have a period of 
3.51±0.40 h. These data are sensitive to puls e fraction^J 
p=U±2% (using eq. 3 from IStute et al.ll201ll l. In Figure 
[2J1 we show the pn light curve in the 0.3-10.0 keV energy 
range phased at the 3.51 h period. We choose the ephemeris 
To=HJD 2454514.54882, which is the date of the firs t pho- 
tometric observation reported bv lWoudt et al.l ((20091). Due 
to the period uncertainty, our phase scale is arbitrary and 
cannot be compared with previous works. 

The U, B and V light curves were scanned for period- 
icities ranging from 1 s to 5 hours. Fourier transform and 



^ We define the pulsed fraction p as max/min = (A+B)/{A—B), 
where A and B are determined from a sine wave fit of the form 
A + B X sm[27r((7i3.5i — </>o)], where 03. 51 is the phase at the 3.51 
h period. 
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Figure 1. Power spectrum Zf{f) from the source-|-background 
(black) and background (red) photon arrival times from the 
pn-|-MOSl,2 cameras. The 99.73% (Scr) confidence level is marked 
by the horizontal line (red). Vertical dotted, dash and dash-dotted 
lines show the 3.94, 3.78 and 3.51 h respectively. A color version 
of this figure is available online. 



phase dispersion minimization algorithms were used in or- 
der to detect sinusoidal as well as complex pulse profiles. No 
pulsed signal with semi-amplitude above 0.03 mag could be 
found in the 1 to 1000s period range in U, B and V bands. In 
particular, no significant peak appears at the 57 s period at- 
tributed to the white dwarf spin; howev er, we note that the 
pulse amplitude of 4 mmag measured by ijWoudt et al.ll2009n 
is below our detection limits. The most significant structures 
in the mid-frequency range are located in the range between 
27 and 32 minutes. Their coherence could not be verified 
and they may just represent the flickering time-scale in this 
system. After rebinning to 100 s integrations, a 4.0±0.8 h 
period is clearly seen in our light curves, with one complete 
cycle covered in U, B and V on different occasions. Simu- 
lations show that a sinusoidal semi-amplitude of 0.13, 0.10 
and 0.14 mag can be measured in U, B and V bands, re- 
spectively, for a 4.0 h hour trial period. Phase-folded light 
curves in the OM filters shows significant departures from a 
sinusoidal shape that may be due to flickering activity dur- 
ing the cycle sampled by our observations (Figure [2]). The 
light curves in each filter show significant change between 
cycles at particular phases. These changes appear as regions 
of larger scatter in the phase-binned curves. 

The individual OM light curves sample less than two 
photometric cycles. Therefore, the broad peak at low fre- 
quencies (with maximum bet ween 3 and 4 hours ) includes 
the 2.889 h period claimed bv lWoudt et al.l l|2009l ) at a high 
power level, producing reasonable phase-folded light curves. 
However, longer optical datasets are needed to study the 
orbital/superhump modulations and the low frequency do- 



3.2 Spectrum 

The X-ray spectrum of V842 Cen is relatively hard with 
photon energies up to > 10 keV (Figure[3| . Visual inspection 
shows an emission feature at energies ~ 6.7-7.0 keV which 
corresponds to Fe XXV and Fe XXVI. A simple thermal 
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Figure 3. XMM EPIC background-subtracted spectra (pn in 
black, MOS 1 in red and MOS 2 in green). The best-fit absorbed 
cooling flow model is overplotted in the top panel. Bottom panel 
shows fit residuals. A color version of this figure is available on- 
line. 



model (with the switch parameter equal 2, meaning that 
the model spectrum used the APEC emissivities) modified by 
a simple absorber (wabs). The model fit (x =260, 267 d.o.f.) 



results in a maximum plasma temperature kTmax = 43^ 
keV, a minimum temperature kT, 



12 

mtn < 0.6 keV, absorption 
, abundance A ^ 2 



column nu = 0.30lo;o4 x 10^^ cm 

solar units and mass accretion rate M[10~^''1=7 _tJW(7) 



m 

Mo/yr 



(d/1.3 kpc)2 (using the distance from lGill fc 0'Brien|[l99i ') 
At the accretion shock front, the shock temperature is 
proportional to the square of the free-fall velocity, which in 
turn is set by the mass of the white dwarf. Thus the white 
dwarf mass can be derived from the value of kTmax obtained 
from fitting the s pectrum with a cooling flow model (e.g., 
lYuasa et aLll2010f ). The value of kTmax obtained from our 
fit implies Mwd=0-9>8 M©, where Mwd is the mass of the 
white dwarf. 



Figure 2. (a) EPIC/pn light curve in the 0.3-10. koV energy 
range folded at the 3.51 h period. Normalized U(6), B(c) and 
V(d) magnitude light curves folded at the 3.51 h period. A sine 
wave fit of the form A + B X sin[2TT{ip3.5i — </'o)]i where 4>3.5i is 
the phase at the 3.51 h period, is overplotted (dashed line). Each 
point corresponds to an average over a 0.05 phase bin. A color 
version of this figure is available online. 



model (x^=297, 271 d.o.f.) of an absorbed plasma (wabs x 
APEC in XSPEq3) yields a temperature kT = 8.lli;g keV, 
an absorption column density uh ~ 0.25^qq3x10^^ cm~^, 
and abundance A = 2.0lQg in solar units. 

Given the detection of a modulation in the X-ray hght 
curves, we also fit the spectrum with models successfully ap- 
plied to cooling accretion shock columns, in which the tem- 
perat ure distribution depends on the cooling mechanisms 
fe.g.. lAllan et aLlll998l ). In the simplest scenario, the X-ray 
spectrum can bo model ed with a n isobaric cooling flow (e.g., 
iMukai et all [2003.: Paridel et ahllioOS. ). We use the mkcflow 



^ http://heasarc. gsfc.nasa.gov /d ocs/xa nadu/xspec/; 
models from ismith et al.l 1I2OIII') : 



iMorrison fc McCammonI l ll98a) 



wabs model 



APEC 
from 



4 DISCUSSION AND CONCLUSIONS 

The nature of V842 Cen remains unclear. Optical data 
obtained in 200 8 show modulatio ns at 56.825 s and 
3.780±0.004 h (|Woudt et all |2009| ). These authors at- 
tributed the shorter period to the white dwarf spin. The 
presence of sidebands at 56.598 and 57.054 s suggested that 
the orbital period should be at 3.94 h. In this picture, they 
then attributed the 3.780 h period to a "negative super- 
hump," a hump-shaped feature in the light curve with pe- 
riod a few percent less than the orbital period of the binary. 
Superhumps are believed to be associated with a tilted ac- 
cretion disk. This analysis strongly suggests that V842 Cen 
should be classified as an IP. It is somewhat puzzling why 
earlier fast optical photometry of V842 Cen did not show 
periodicity (jWoudt fc W arner 20031); presumably, flaring ac- 
tivity observed in that data masked any periodic signatures. 
Our data cannot confirm the short period at X-ray or optical 
wavelengths. While the lack of a ~57 s period in the X-ray 
data does not necessarily rule out a rapidly rotating white 
dwarf, one generally expects the X-ray data from IPs to be 
modulated at the spin period, since the localized accretion 
shock near the surface of the white dwarf produces a high 
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temperature thermal spectrum. Unfortunately, our optical 
data are not sensitive enough to detect a pulse amplitude of 
4 mmag as previously measured. 

On the other hand, we detect significant periods at 
3.5f ±0.40 h in X-rays and 4.0±0.8 h in optical. The pulsed 
fraction p of the X-ray light curves folded at this period de- 
creases with energy, ranging from 47±13% in the soft (0.3- 
1.0 keV), to 45±12% in the medium (1.0-2.0 keV), and to 
30±11% in the hard (2.0-10.0 keV) X-ray band. A similar 
trend is observed in the X-ray spin-phased light curves of IPs 
(e.g.. lAllan et al.lll998l . EX Hya) and can be caused by the 
energy dependence of the absorption cross-section, occulta- 
tion of the lower (cooler) portion of the accretion column by 
the white dwarf's body or both. The modulation detected 
in our data is somehow different from the orbital period 
claimed bv lWoudt et al.l (|2009l ) at the 3(t level. Furthermore, 
the smooth shape of the phase-folded light curves indicates 
that an orbital origin of the mod ulation is unlikely. A low 
orbital inclination as suggested bv lWoudt fc Warneij (12003) 
also supports this conclusion. Eclipses of the accretion col- 
umn are observed as sharp decrea ses in the light curves o f 
IPs with high orbital inclination iJHoogerwerf et al.ll200a ). 
lasting for tenths or less of the orbital period. No conspic- 
uous eclipse feature could be found using a 50 s binning in 
any of the light curves. 

On the other hand, the period we have detected is con- 
sistent with the value reported bv lWoudt et al.l (|2009l ) as a 
negative superhump; however, it is not clear how the tilted 
accretion disk would contribute to such high energy X-ray 
emission. To our knowledge, negative superhumps have not 
been detected in a CV at X-ray wavelengths. Furthermore, 
the energy dependence of the pulsed fraction of an X-ray 
superhump would require explanation. 

The X-ray spectrum is fully consistent with thermal 
emission originating in the post-shock regi on of a magnetic 
accretion column (e.g. iMukai et al.l 120031 ). The model fit 
to the X-ray spectrum suggests that the white dwarf in 
V842 Cen has a mass Mwd=Q-?>^ Mp,, whi c h is n ot un- 
usual for an IP. In fact, iBrunschweiger et al.l (|2009l ') found 
that among the IPs detected with Swift/HAT, 5 out of 
17 have white dwarfs more massive than 0.88 Mq. Given 
the good agreement with the accretion shock model and a 
reasonable mass determination, it is difficult to reconcile the 
lack of detection of the fast period associated with the white 
dwarf spin. Either the fast period is undetectable due to a 
low inclination or the period we detect is instead the white 
dwarf spin period. The latest possibility would be, if con- 
firmed, very rare since the detected period is longer than 
any known spin periods from IPsq 

The shape of the nova remnant in V842 Cen also re- 
mains puzzling. The remnant geometry depends on the mat- 
ter distribution and the ejecta illumination by the central 
source. If no accretion disk is present in V842 Cen (e.g. it 
is a polar CV instead of an IP), then the ejecta might be 
more homogeneously illuminated; however, at the moment 
we do not think that this is a possibility as no evidence has 
been reported elsewhere on the detection of optical polar- 
ization or strong soft X-ray emission that would support the 



^ See the latest catalog of IPs at 

|http://asd.gsfc.nasa.gov/Koji.Mukai/iphome/catalog/alpha.html| 



polar nature. One expects the presence of a large accretion 
disk to affect the photoionization in different regions of the 
nebula. The accretion disk may absorb the ionizing photons 
in the equatorial regions causing an aspherical illumination. 
In HR Del for example, the observed aspherical illumina- 
tion is att ributed to high mass transfer and a large disk 
(JMoraes fc Diaz 2009 ). 

Furthermore, a high white dwarf rotational velocity can 
generate asymmetries in the eruption conditions because of 
the eff ective gravi ty variation from equatorial to polar re- 
gions (jScottllioOol ). If we consider a 0.88 Mq white dwarf 
with a rotational period of 57 s, we should obtai n ejecta wit h 
axial ratio of ~1.35 (using expressions 8 to 11 bv lScottl200d '). 
If the rotational period is 3.51 h, the envelope would be sym- 
metric. A period P<125 s is required to detect an asymmet- 
ric envelope. If the white dwarf is more massive than what 
we have derived, the asymmetry would be smaller, but de- 
tectable (as a prolate remnant). 

In summary, the X-ray spectrum and the energy de- 
pendence of its pulse fraction support the classification of 
V842 Cen as an IP, but the symmetric shape of the ejecta 
and the long period of the light curves are not easily ex- 
plained in this picture. While the X-ray data do not defini- 
tively rule out a very fast 57 s spin period, a longer period 
would be less discrepant with the ejecta symmetry. Another 
possibility is that the ejecta asymmetry axis is the same as 
that of the line of sight. A clear, definitive determination of 
the nature of V842 Cen, in particular whether or not it is 
an IP and what its spin period is, will contribute to under- 
standing the system and its ejecta morphology. 
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